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Diquark excitations of dense quark matter are considered in the framework of the Nambu - 
Jona-Lasinio model with three types of massless quarks in the presense of a quark number chemical 
potential /j,. We investigate the effective action of meson- and diquark fields at sufficiently high values 
of fi, where the color-flavor locked (CFL) phase is realized, and prove the existence of NG-bosons 
in the sector of pseudoscalar diquarks. In the sector of scalar diquarks an additional NG-boson is 
found, corresponding to the spontaneous breaking of the U(l)s baryon symmetry in the CFL phase. 
Finally, the existence of massive scalar and pseudoscalar diquark excitations is demonstrated. 
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C i I. INTRODUCTION 

d . 

. It is well-known that at sufficiently high baryon densities massless three-flavor QCD is in the so-called color - 
flavor locked (CFL) phase 0,0. In this phase the original SU(3) L xSU(3) fl xSU(3) c xU(l) s symmetry of QCD is 
spontaneously broken down to the diagonal SU(3)l+_r+ c subgroup. Correspondingly, seventeen massless excitations 
£SJ ■ must appear in the mass spectrum of the theory. Eight of them might be used to ensure a mass of gluons by the 
^ [ Anderson - Higgs mechanism. The properties of the remaining nine, one scalar and eight pseudoscalar Nambu - 
• Goldstone (NG) bosons.as well as other collective modes of the CFL phase were studied already in the framework of 
OO ; weak-coupling QCD @, i, [| . 

It is clear that a weak-coupling QCD analysis of the color superconductivity phenomena, including the CFL one, 
can only be trusted at rather high baryon densities or, correspondingly, for values of the quark number chemical 
potential /j, >> 1 GeV [6]. At moderate values of fx ~ 500 MeV, where weak-coupling QCD is not applicable, the 
description of color superconductivity is more adequate in the framework of effective theories for the low energy QCD 
region. In particular, since massless excitations might play an important role in different transport phenomena, a 
chiral effective meson theory for the pseudoscalar NG-bosons of the CFL phase was also constructed 0, H| (see also 
I ' the recent review on this topic [9[ and references therein). 

Another effective quark model approach to low density QCD is based on the Nambu - Jona-Lasinio (NJL) type of 
models. Since the NJL model contains the microscopic quark degrees of freedom, it is more fundamental and preferable, 
especially for the investigation of dynamical processes in dense baryonic matter, than chiral meson theories. On the 
other hand, it is also more suitable for the description of physics at low baryon densities than weak-coupling QCD. 
In p articular, in the three-flavor NJL model the CFL effect was already considered, e.g., in (see also the review 
[111]), where some aspects of the phase structure of dense quark matter were discussed, including the influence of the 
s-quark bare mass, color- and electric charge neutrality conditions, external magnetic field, etc. However, in spite 
of the fact that the lightest bosons may play an essential role e.g. in the cooling processes of neutron stars, up to 
now only few attention was paid to the consideration of the CFL ground state bosonic excitations, i.e. mesons and 
diquarks, and their dynamics in the framework of NJL models (see, however, the papers [l2l.[l~3l|. where the properties 
of the massless NG boson, corresponding to the spontaneous breaking of the baryon U(l)s symmetry in the CFL 
phase, were considered). In contrast, the properties of w- mesons and diquarks, surrounded by color sup erconducting 
quark matter, were already discussed in the framework of the two-flavor NJL model 0, [H, EE US, ITa | ■ 

In the present paper we are going to study meson and diquark excitations of the CFL ground state in the framework 
of the massless three-flavor NJL model. As in [3, EH EH, our consideration is based on the effective action, which 
is a generating functional for one-particle irreducible Green functions. They permit to get informations about the 
masses of bosonic excitations of the CFL phase. In the first step we will rederive the well-known result that in the 
normal phase of quark matter the NG bosons are just pseudoscalar mesons. Next, it will be shown that in the CFL 
phase the NG bosons are scalar- and pseudoscalar diquark excitations. Finally, we will demonstrate that in the CFL 
phase there appear massive scalar and pseudoscalar diquarks which are composed into a charged triplet and singlet 
as well as a neutral singlet of the SU(3) group. 
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II. NJL MODEL AND EFFECTIVE MESON-DIQUARK ACTION 

Let us consider the following NJL model with three massless quark flavors 

8 

L = qfridv + A*7° q + G^ [(qr a q) 2 + {qil 5 r a q) 2 + 

a=0 

+ G 2 ^2 ^ [[q C ij 5 T A \ A ,q][qij 5 T A \ Al q C ) + [q C T A \ Al q][qT A \ Al q C ]y 
A=2.5,7 A'=2,5,7 



(1) 



In ([T]), /i > is the quark number chemical potential, which is the same for all quark flavors, q = C(f, q c — q t C are 
charge-conjugated spinors, and C — ij 2 ^f° is the charge conjugation matrix (the symbol t denotes the transposition 
operation). The quark field q = qi a is a flavor and color triplet as well as a four-component Dirac spinor, where 
i,a = 1,2,3. (Roman and Greek indices refer to flavor and color indices, respectively; spinor indices are omitted.) 
Furthermore, we use the notations r Q ,A a for Gell-Mann matrices in the flavor and color space, respectively (a = 

1, ...,8); To = 1 is proportional to the unit matrix in the flavor space. Clearly, the Lagrangian ([T]) as a whole 

is invariant under transformations from the color group SU(3) C . In addition, it is symmetric under the chiral group 
SU(3)lxSU(3).r (chiral transformations act on the flavor indices of quark fields only) as well as under the baryon- 
number conservation group U(l)s and the axial group U(1)a- 1 

The linearized version of the Lagrangian ([I]) contains collective bosonic fields and looks like 



L 



q 

A 



j"id„ + ^7° - a a r a - ij 5 Tr a T a 



1 



1 

AG~2 



_1_ \P* AP 
A A' + AA' AA' 



AA 



-[q C i>y 5 T A \ A/ q} 



A 



AA 



-[qij 5 T A X A >q c 



AP* 

^[q C T A \ A ,q] 



^[qr A \ A ,q% 



(2) 



where here as well as in the following the summation over repeated indices a = 0, 8 and A, A' = 2,5,7 is implied. 
Lagrangians ([T]) and are equivalent which simply follows from the equations of motion for the bosonic fields 



a a (x) = -2Gi(qT a q), A s AA ,{x) = -2G 2 (q c i~f 5 T A \ A ,q), A s A * A ,{x) = -2G 2 (^ 7 5 r A A^( Z G ), 
n a (x) = -2G 1 (^ 7 5 r a ( 7 ), A AA ,(x) = -2G 2 (q C T A X A ,q), A p A * A ,(x) = -2G 2 (qT A X A ,q c ). 



(3) 



One can easily see from §S§ that the mesonic fields a a (x) , ir a (x) are real quantities, i. e. (a a (x))' i = (T a (x), (7r a (x))t = 
Tr a (x) (the superscript symbol f denotes the hermitian conjugation), whereas all diquark fields A^, (a;) are complex 
ones, i. e. 



MA' 



,(x), (A* AA ,(x))<< = A A * A ,(x) 



Moreover, A s AA ,(x) and A AA ,(x) are scalars and pseudoscalars, correspondingly. 

Let us introduce the flavor group SU(3)/=SU(3)l+^, which is the diagonal subgroup of the chiral group. Then, 
all scalar diquarks A s AA ,(x) form an (3 C , 3/)-multiplct of the SU(3) c xSU(3)/ group, i. e. they are a color and flavor 
antitriplct. The same is true for pseudoscalar diquarks A AA ,(x) which are also the components of an (3 C , 3/)-multiplct 
of the SU(3) c xSU(3)/ group. Evidently, all diquarks A AA ,(x) have the same nonzero baryon charge. All the real 
a a (x) and ir a (x) fields are color singlets. Moreover, the set of scalar cr a (a;)-mesons is decomposed into a direct sum of 
the singlet and octet representations of the diagonal flavor group SU(3) t . The same decomposition into multiplets is 
true for the set of all pseudoscalar 7r a (a;)-mesons. Clearly, in this case the octet is constructed from three pions (it 
and 7T°), four kaons {K°, K° and K^) and the eta-meson (?y 8 ), whereas the singlet (770) corresponds to the 77 '-meson. 

Next, in order to use the Nambu-Gorkov formalism, we put the quark fields and their charge conjugates together 
into a bispinor ^ = ( c) so that the Lagrangian ([2|) takes the compact form 



L = — 



4Gi 



1 

4G^ 



2 



where Z is the 2 x 2-matrix in the Nambu-Gorkov space, 



-K* 



-K 
D- 



(4) 



(5) 



1 In a more realistic case, the additional 't Hooft six-quark interaction term should be taken into account in order to break the axial 
U(1)a symmetry Il0( . However, in the present consideration we omit the 't Hooft term, for simplicity. 
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and the following notations are adopted 



D+ = i-fd v + ^7° - S, E = r a a a + ij 5 ir a T a , K = {A P AA , + iA AA ,~/ 5 )T A \A> , 
D~ = i^ft, - ^ - E«, E t =r* ( r a + z7 5 7r Q T*, K* = (A^, + iA aa ,j 5 )taXa' 



(6) 



Now, integrating out the quark fields in the partition function with the Lagrangian |4]) (this procedure is presented 
in detail in [l6[), we obtain the effective action for the considered NJL model in the one-fermion loop approximation, 



dS 



tt; A s aa ,A aa , + A AA ,A A * A , 



4Gi 



+ 



4G 2 



Tr In Z. 

2 



(7) 



Besides of an evident trace over the two-dimensional Nambu-Gorkov matrix, the Tr-operation in ([7]) stands for 
calculating the trace in spinor, flavor, color as well as four-dimensional coordinate spaces, correspondingly. Let us 
suppose that parity is conserved so that all pseudoscalar diquark and meson fields have zero ground state expectation 
values, i. e. (A^^a;)) = and ('K a {x)) = 0. Furthermore, since at zero s-quark mass, m s — 0, only the competition 
between the normal quark matter phase and the CFL one is relevant to the physics of dense QCD (see, e.g., Fig. 1 in 
0, where the corresponding phase portrait for QCD with two massless u, <i-quarks at zero temperature is presented 
in terms of /i, m s ), we permit in the present consideration nonzero ground state expectation values only for ao(x) and 
A s AA (x) fields (A = 2,5,7). Namely, let (a (x)) = erg, (A s AA (x)) = A, (A s A * A (x)) = A*, where A = 2,5,7, but other 
boson fields from ([3]) have zero ground state expectation values. In this case, if A = 0, then quark matter is in the 
normal phase, where at erg ^ the ground state is an SU(3) c xSU(3)/xU(l)_B-invariant one. If A ^ 0, then the CFL 
phase is realized in the model, and the initial symmetry is spontaneously broken down to SU(3)l+^+ c . Now, let us 



make the following shifts of bosonic fields in ([7]): er (a 



<7 (x)- 



A A * A (x)+A*, A AA (x) - A AA (x)+A, 



where A = 2, 5, 7, and other bosonic fields remain unshifted. (Obviously, the new shifted bosonic fields ao(x), A s AA (x) 
etc, now denote the (small) quantum fluctuations around the mean values erg, A etc of mesons and diquarks rather 
than the original fields ©.) In this case 



D+ , -K 
~K* , D- 



S , K\ =1 
K* , S* : " 



E 

K* 



K 
E* 



(8) 



where K a , K*, Df, E Q , E* are the corresponding quantities ([6]), m which all bosonic fields are replaced by their own 
ground state expectation values, i.e. oq(x) — > erg, ir a (x) — > 0, A s AA (x) —> A, A P AA , (x) — ► etc, and So is the quark 
propagator matrix in the Nambu-Gorkov representation (its matrix elements are given in the Appendix) . Then, 
expanding the obtained expression into a Taylor-series up to second order of small bosonic fluctuations, we have 



Seff (cr a , 7r a , A S AA , , A^ ) = S$ + S ( J> (a a ,ir a ,A" AA ,, A s / A * ) 



:(2) 



(9) 



where 



c(0) _ 



S$(a a ,n a ,A s AA ,,A s AA *) = - I d*x 



4Gi 



3|A| 

4G 2 



|Trln(S ) = -nK,A,A*) J d 4 



+ 4 Tr i 5 ° { K*\ E* 



a 2 „+n 2 „ A AA ,A A \,+A p AAI A p ; A , 



AA' AA' 



4Gi 
So 



4G 2 



E , K 
K* , E* 



+ 



(10) 



(11) 



Notice that the term linear in meson and diquark fields vanishes in © due to the gap equations (see below). In |T 
the quantity f2 is the thermodynamic potential of the system. In terms of M = \i |erg it looks like 



n(M, A, A*) 



3M 2 



3AI* 



8Gi 4G 2 



d 3 q 
(2tt) 3 



E 2A 



E. 



2A 



Y 



(12) 



where (E^) 2 = (E*) 2 + |A| 2 , {E± A ) 2 = {E^f + 4|A| 2 , E ± = E ± fi, E = J q 2 + M 2 . Starting from ([12]), one can 

find the gap equations dil/dA* = and dQ,/ dM — 0, which supply us with the values of M, A in the ground state 
of the system: 



dM 

dn 

dA* 



3M 
4GT 
3A 

: 4G 2 



2M 



d 3 q f 4E+ 
d 3 q f 4 4 

j^r\Ej + E^ 



4E- 



E+ 

2E 2A 



E. 



E' 



2E. 



2A 



o, 



2A 



(13) 
(14) 



4 



M 




FIG. 1: The behaviour of the thermodynamic potential Q. 
vs M and A at the critical value of the chemical potential 
H C « 329 MeV. 



FIG. 2: The coordinates M and A (gaps) of the global 
minimum point of the thermodynamic potential vs the 
chemical potential /i. Here fi c w 329 MeV. 



Since the integrals in the right hand sides of l|12 |) -(|14 p are ultraviolet divergent, wc regularize them as well as the 
other three-dimensional divergent integrals below by implementing a cutoff in the integration regions, \q\ < A. In all 
numerical calculations below, we use the following values of the model parameters (see, e.g., ref. [ill ]) 

A = 602.3 MeV, dA 2 = 2.319, G 2 = 3<?i/4. (15) 

Suppose that A is a real quantity. In this case, since the thermodynamic potential (fT2"|) is an even function in 
both the M and A variables, it is enough to study it only in the region {M > 0, A > 0}. Just for this region 
a graph of (l(M, A) is presented in Fig. 1 at fi = fx c « 329 MeV. It is clear from Fig. 1 that in this case the 
thermodynamic potential has two global minimum points (GMP), A and B. The point A with coordinates (M c , 0), 
where M c w 355 MeV, lies on the M-axis and corresponds to the normal quark matter phase, whereas the point 
B with coordinates (0, A c ), A c « 86 MeV, is arranged on the A-axis and corresponds to the CFL phase of dense 
baryonic matter. (Evidently, both A and B are solutions of the gap equations (fT3]) -(114 |) . In addition, as is seen from 
Fig. 1, there are another solutions of the gap equations which, however, correspond to a maximum- or saddle points of 
the thermodynamic potential.) If /i < /i c , the function fi has only one GMP of the A-type, i.e. it lies on the M-axis. 
So, in this case the normal quark matter phase is realized in the model. In contrast, if \i > /z c , then the single GMP 
of the fi-function is located on the A-axis, i.e. it is of the -B-type, and the CFL phase occurs. The dependence of the 
coordinates (gaps) of the GMP on the chemical potential is presented in Fig. 2. Since the GMP jumps at /i = /i c 
from A to B (or vice versa), one may conclude that at /i = /i c a first order phase transition takes place in the system. 

The effective action S^g in is really a generating functional of the one-particle irreducible (1PI) two-point 
Green functions of mesons and diquarks both in the normal and CFL phases, namely 

r ^ x -y^-im^y (16) 

where X(x),Y(x) — a a (x), 7Tb(x), A^, (x), A^*, (x). These Green functions are very useful, in particular, in de- 
termining the dispersion relations and masses of particles. In the following, we shall say that in the theory there 
is a mixing between two different particles with corresponding fields X(x) and Y(x), if their 1PI Green function 
Lxy(x — y) is not equal to zero. By analyzing the structure of the effective action (jlip . it is possible to show that 
in the NJL model ([1]) with three massless quarks, mesons and diquarks are not mixed both in the normal (fx < [i c ) 
and CFL phases {fi > fi c ). 2 Moreover, each pseudoscalar 7r (x) field (as well as scalar <7&(x) one) does not mix with 
other meson fields, apart from itself. So, in the normal phase (where A = and M ~ 355 MeV) one can obtain 
the following expression r,^,^ (p) for the Fourier transformation of the two-point 1PI Green function of pseudoscalar 
mesons Wi(x), taken in the rest frame p = (po, 0, 0, 0): 

. . 2 f d 3 q 12 
TwM=PoJ E[p 2_ AE 2 Y (!7) 



2 Note, if some of the current quark masses are nonzero, then there arises a mixing between mesons and diquarks in the CFL phase. This 
effect is the analogy of the mixing between the cr-meson and the scalar diquark in the color superconducting phase of a two-flavor NJL 
model with nonzero masses of u- and d-quarks Il5l . I 111 . 
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(In obtaining (| 1 T[) . the gap equation (fT3|) was used in order to eliminate the coupling constant G\ from the final 
expression.) Evidently, this expression turns into zero at p\ = 0. Since the relation (| 1 7[) is true for i = 0, 1, ...,8, it 
means that nine massless excitations, Nambu - Goldstone (NG) bosons, do exist in the pseudoscalar meson sector 
of the model in the normal phase. This fact corresponds to a spontaneous symmetry breaking in the normal phase 
down to the group SU(3) c xSU(3)/xU(l) s . 

Recall, that our initial NJL model (p} is invariant under the symmetry group SU(3)l xSU(3)_rxSU(3) c x 
U(1)b xU(1)a, which contains an additional (nonphysical) axial subgroup U(1)a- So in the CFL phase, where the 
ground state is symmetric under the group SU(3)l+_r+ c , an additional NG boson must appear in the mass spectrum 
of the model giving altogether eighteen NG bosons. 

Note that in the sector of the meson excitations of the CFL phase we do not find any massless particle. 



A further consideration of the effective action (fTTj) shows that in the CFL phase, where M = and A ^ 0, 
scalar and pseudoscalar diquarks are separated from each other. Moreover, the eighteen pseudoscalar diquarks (nine 
A AA ,(x) an d nme &?AA'( X ) ne lds) may be divided into four sectors: p(57, 75), p(25, 52), p(27, 72) and p(257). Each of 
the sectors p(AA',A'A), where A ^ A\ is composed of A P AA , (x), A A * A , (x), A p A , A {x) and A^* /j4 (x)-diquarks, whereas 
the sector p(257) is composed of six fields, A^ix), A P l(x), A P j(x), A p 2 (x), A p 5 (x), and A p 7 (x). There is a mixing 
between diquarks from the same sector, however fields from different sectors are not mixed. (The analogous situation 
takes place for the set of scalar diquarks.) 

Partially, these mixing properties of the pseudoscalar diquarks are explained by the ground state SU(3)l+t?+ c 
symmetry of the CFL phase. Indeed, with respect to this group all pseudoscalar diquarks are decomposed into a 
direct sum of the 3 and 6 multiplets. The fields A p AA (x) with A = 2, 5, 7 are in an antitriplet 3, whereas all diquarks 
of the form A p A , A {x) (A 1 ^ A) are in an antisixtet 6. It is clear from symmetry considerations that diquarks from 3 
do not mix with diquarks from 6 (here we use the terminology, introduced just after (|16|) ). Since the mixing between 
a particle and corresponding antiparticle is allowed as a rule, we may conclude that pseudoscalar diquarks from the 
sector p(257) are separated from the other six A A , A {x)- and six A^*, A (x)-diquarks [A' ^ A). A further separation 
between components of the 6-multiplet and corresponding antiparticles occurs on a dynamical basis, i. e. it is due 
to the structure of the effective action (jlip and the quark propagator matrix (|A.ljl - (|A.4j) . As a result, one can show 
that these diquarks are divided into the three above-mentioned sectors, p(57, 75), p(25, 52), p(27, 72). 

Below we suppose that the gap A is a real positive number in the CFL phase. 



Let us first study the mass spectrum of the excitations, e.g., in the sector p(57, 75). The two-point 1PI Green 
functions of pseudoscalar diquarks from this sector can be obtained from by the relation (fT6| . In the rest frame, 
where p = (po, 0, 0, 0), the Fourier transforms of these 1PI Green functions form the following matrix: 



III. MASSLESS DIQUARK EXCITATIONS IN THE CFL PHASE 



A. The case of p(AA' , A' A) sectors 




After tedious but straight-forward calculations, based on the technique elaborated in our previous papers 0, EH, EH 
(see also [l7], [U|) an( i used in the consideration of mesons and diquarks in the color superconducting phase of the two 
flavor NJL model, the matrix (fT8]l takes the form: 




(19) 
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where A = a + pop, B = a — pa[3 and 

d 3 Q (6E+p 2 -(E+ + E+ A f(2E++E+ A ) , 4p 2 + A(E+ A ) 2 - 10A 2 



P 



(2tt)3 1 9EXE+ A [pl (E+ + E+J*] SE+lpl 4(25+)»] 

d 3 q 
(2tt) 3 

d 3 g f £+ A ) 10E^ 



~* E A J ^2A ~* r ' ( 20 ) 



(2tt)3 1 3£+i? 2 + A b 2 - (£+ + £+ A ) 2 ] 3£+b 2 - 4(£+) 2 ] 



rf 3 g 
(2tt) 3 



, E A ^>E A , E 2A ^E 2A \, (21 



r , d\ f 2A 2 (£+ + g+ A ) 10A 2 

(2^) 3 \ 3£+£+ A b 2 - (J5+ + £+ A ) 2 ] 3£+b 2 - 4(£+) 2 ] 



+ j J^nf ^sA^^Aj- ( 22 ) 

(To obtain the above expressions for a and /3, we have used the gap equation (|14p in order to eliminate the coupling 
constant Gi from corresponding 1PI Green functions.) Evidently, in the case of mixing between particles the infor- 
mation about their masses should be extracted from zeros of a matrix determinant, composed from corresponding 
1PI Green functions. So, in our case it is necessary to study the equation detl^^^o) = 0, which takes the following 
form 

detr 57 ,7 5 (Po) = (AB - C 2 ) 2 = [(a - C)(a + C) - p 2 p 2 } 2 = 0. (23) 

In the pp-plane, each zero of this equation defines a mass squared of a bosonic excitation of the CFL phase ground 
state in the p(57, 75) sector. Since this sector contains four pseudoscalar diquarks, one should search for four solutions 
of the equation (|2"3"|) in the p^-plane. However, due to the structure of (|23|) . it is clear that this equation admits at 
least two different solutions (each being two-fold degenerate) , which are given by the zeros of the expression in the 
square bracket. Since 

„ r-^ f d3q fPll 2E A+^+ A }-8(E+f-(E+) 2 E+ A -10E+(E+ A ) 2 -5(E+ A f 



(2^) 3 I 3E+E+ A [p 2 (E+ + E+ A n P 2 4(£+) 2 ] 

+ Plj 7^X3 {^a - E A , E+ A - E- A | = plF{pH (24) 

the square bracket in ([23]) becomes zero at the point p 2 , = 0. So, in the p(57, 75) sector there are two massless 
pseudoscalar excitations, i.e. NG bosons. Two other excitations have the same nontrivial mass squared which is the 
solution of the equation 

(a + C)F(p 2 )-p 2 = 0. (25) 

Its investigation is outside the scope of the present paper. A similar situation occurs in the other four-component 
sectors p(2b, 52) and p(27, 72). Namely, for both sectors the 1PI Green function matrix has the form (fT9|) . Hence, in 
each of these sectors there are two NG bosons as well as two massive excitations with the same mass squared. Its 
value is given by the solution of the equation (|25l) . 

To summarize, in the pseudoscalar diquark sectors p(57, 75), p(25,52), and p(27, 72) we have found six massless 
diquark excitations, which are NG bosons, as well as six pseudoscalar diquark excitations with a nonzero mass which 
is a solution of the equation (|25|) . In total, these six massive real pseudoscalar diquarks form a complex (charged) 
triplet of the SU(3)l+,r +c group. 



The case of the p(257) sector 



All the two-point 1PI Green functions of the pseudoscalar diquarks, entering the p(257) sector, are defined by the 
relation (|16p . In momentum space representation and in the rest frame, i.e. at p = (po, 0, 0, 0), they form the following 
6x6 matrix: 



^57 (?<)) = 



r 
r 

—55 — 

r A^A?,(Po) 



AS>o) 

Af 5 Af 2 (Po) 



_A? r A| 



V r A? 



>o) 
>o) 



r A f 2 A§ 2 *(po) 

£Af 5 Af 2 *(Po) 

£as;a;;(po) 

£a? 7 a-(po) 

r A? 7 *A^(Po) 



£Af 2 Af 5 (Po) 

£ A | 2 * A f 5 (po) 

£a* 5 a{; 5 (Po) 
£aI 5 *a| 5 Oo) 

£A? 7 Ag 5 (Po) 



i»;(po) 

A* 5 Ag*(Po) 
A|*A|*(Po) 

^ 7 a?:(po) 

77 55 



Af 2 A? 7 (Po) 

Ag 5 A? 7 (Po) 
?JA- 7 (P0) 
7 A? 7 (Po) 



r 



As, A 



a?: a 



f B (po) r A?;A g 5 .(po) 



a; 

i_i 77 AA. 



»(Po)\ 

r A f*A? 7 *bo) 

£a? 5 a?;(po) 

£aj;a?;(po) 
r 



A?,A 



::Cpo) 



r (po) r A? 7 *A? 7 *bo) / 



(26) 
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A tedious but straight-forward calculation yields: 



T257(po) = 



/ R P T W T W\ 

Q R Z T Z T 

T W R P T W 

Z T Q R Z T 

T W T W R P 

\ztztqrJ 



The determinant of this matrix looks like 

detr 257 (po) = [(T - Rf -{W- P)(Z - Q)} 2 {{2T + R) 2 - (2W + P)(2Z + Q)}, 



where P = I + p h , Q = h - Poh , W = J 

d 3 q 



PoJi, Z = J -p Ji, and 



h(p 2 ) 



28 

E ni ^)[E A E 2A 



(2tt) 3 
+ 

2A 



1 

2 (El 

9 E+E+ A [p 2 -(E 

d 3 q f 28 



A) 



(E+f 



4 (E+ A ) 2 + (E+) 2 
EX\pl-4(E+r] ' 9 E+ A {pl-4(E+ A )2] 
{E+) 2 ]} , f d 3 q 



E+ 



(27T) 



E A,2A 

75+ 



^A,2A' 



(2^) 3 l9 4(^) 2 ] 9 £+ A [p 2 -4(£+ A ) 2 ] 

(E++E+ A )E+ 1 , /■ A 



9 EZE+^tf 



R( P 2 ) = (— A ) 



A 2A U'O 

J 3 q 



A ^ 

28 



£ 2 + a) 2 ] 



(271)3^9 S+[pg-4(S+)2] 

d 3 



(2tt) 3 1 A < 2A 
16 



A,2A 
1 



E+ A [p 2 ~A(E 2A 



E+ + E. 



2A 



•/o(pg) 



9 EiE} A [p 2 ^(E++E+ A ) 2 } 

d 3 q (4 _(E+ A ) 2 + (E+) 2 
9 



9 



(2tt) 3 
2 (£ 



A,2A 



i) 2 



2ALP0 
— * ^A,2A 

(E+f 



(2tt) 3 19 E+ A [p 2 -A(E 2A 

i ( J B+ + i? 2 + A )[£;+i?+ A + (ii;+) 2 ] 

9 



E A E 2A 



Jlipl) 



d 3 q 



4 



^2 + A) 2 ] 



+ 



9 i&[p§-4(i5£)*] 
J 3 q 



(2tt) 3 



^A,2A 



E 



A,2A 



(2tt) 3 19 £+ A [p 2 -4(£+ A ) 2 ] 9 E+[p 2 -A(Ei) 2 } 
(E++E+ A )E+ \ , /■ A 



1 



9 £+£ 2 + A b 2 



^2 + a) 2 ] 



(2tt) 3 



A.2A 



E 



A,2A 



E~ 



-E~ 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



^) = (-A 2 )/|f,{ 



16 

IT 



^2A 



^2 + a[P0 



E+E+ A \p 2 



J A ±J 2A\-I J ( E A 

In terms of the J;-functions, the expression 



^2 + A) 2 ] 



1 



4(S 2 + A ) 2 ] 
J 3 g 



<2AJ J , - (2^) : 
can be rewritten as 



9 E+[p 2 -A(E+y} 



^A,2A 



E 



A,2A 



(34) 



detr 257 (po) = [(T - i? - Jo + /o)(T -R+J -I ) + p 2 (Ji - h?\ {{2T + R 
- 2 Jo - I )(2T + R + 2 Jo + J ) + pl(2Ji +h) 2 }. 



(35) 



As in the previous section, the full spectrum of the CFL phase excitations in the p(257)-sector of the model can be 
explicitly computed by evaluating the zeros of the determinant ([28)1 (or (f35"j) . alternatively). Since the p(257)-sector 
consists of six diquark degrees of freedom, we expect that defies? (po) has at least six zeros in the Po~plane. 

Eliminating in (|29p the coupling constant G 2 again with the help of the gap equation (|14p , one can easily show 
that T(0) = Jo(0) and R(0) — Iq(0)- Evidently, in this case both the square bracket expression and the brace one 
from (|35|) turn into zero at p 2 , = 0. Hence, in the p(257)-sector of the model there are three pseudoscalar NG boson 
excitations. Now let us obtain some information about the remaining three excitations, which are massive ones. For 
this purpose, similarly to (|2~4"|) , it is possible to extract in an evident form the factor p\ from the expressions in the 
curly brackets of (|35"|) , i. e. 



T - R - Jo + I = pIHpI), 2T + R - 2 J - I = pM P 2 ). 



(36) 
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Then it is clear from (|35| that two of these excitations have an identical mass which is a solution of the equation 

4>(p 2 )(T -R + J q -I ) + (Ji - h) 2 = 0. (37) 

In total, they form a complex (charged) singlet of the group SU(3)a+l+ c . Finally, there is a further diquark excitation, 
whose mass obeys another equation 

<p(p 2 Q )(2T + R + 2J a + Jo) + (2Ji + h) 2 = 0. (38) 

Evidently, this diquark is also a SU(3)i?+L+ c -singlet. As a result, we have proved that there are nine pseudoscalar 
NG bosons (diquarks) in the mass spectrum of the initial NJL model. Moreover, it contains massive pseudoscalar 
diquarks, composed of a complex triplet and singlet as well as real singlet of the group SU(3)/+ C . 

Since in the chiral limit, where all quarks are massless, the two-point 1PI Green functions of scalar diquarks are 
identical to the corresponding Green functions of pseudoscalar diquark fields, one may conclude that the diquark 
spectrum contains nine scalar massless excitations. In real QCD, eight of them should supply a mass to gluons 
through the Anderson - Higgs mechanism, but the remaining massless excitation is the NG boson corresponding to 
the spontaneous breaking of the baryon U(l)s symmetry. In addition, there are a massive complex triplet and singlet 
as well as a real singlet of scalar diquarks in the mass spectrum of the model. 

IV. CONCLUSIONS 

In the present paper the two-point 1PI Green functions of scalar and pseudoscalar diquarks are investigated 
in the framework of a three flavor NJL model with massless quarks and nonzero chemical potential /i. The 
model contains interaction terms both in the quark-antiquark and quark-quark channels, but the 't Hooft six- 
quark term is omitted, for simplicity (see (H|)). In this case, the initial symmetry group of the model, i. e. 
SU(3)lxSU(3)_r,xSU(3) c xU(1)sxU(1) j 4 does contain the axial U(1)a subgroup. As a result, we have shown that at 
sufficiently low values of /i, < /x c » 330 MeV, the normal quark matter phase with SU(3)l + _rxSU(3) c xU(1)b is 
realized and nine massless pseudoscalar mesons (which are the NG bosons), 7r ± , ir°, K°, K°, K^, 7]$ and 77', appear. 
(In massless QCD, where 11(1)^ is broken on the quantum level, or in NJL models with 't Hooft interaction the 
7] '-meson is not a NG boson.) 

At \x > \i c the original symmetry of the model is spontaneously broken down to SU(3)l+^+ c , and the CFL phase 
does occur. In this case, in accordance with the Goldstone theorem, eighteen NG bosons must appear in the mass 
spectrum. Considering 1PI Green functions, we have found nine NG bosons in the sector of scalar diquark excitations. 
Eight of them have to be considered as non-physical, since in real QCD they supply masses to gluons by the Anderson 
- Higgs mechanism. The remaining scalar NG boson corresponds to a spontaneous breaking of the baryon U(l)s 
symmetry. The other nine NG bosons are no more pseudoscalar mesons, but now the massless excitations in the 
pseudoscalar diquark sector of the model. 

Besides NG diquarks, we have proved the existence of massive diquark excitations in the CFL phase. They form 
both pseudoscalar and scalar complex (charged) triplets and singlets as well as a real (neutral) singlet of the group 
SU(3)/+ C . There arises the interesting question, whether these diquark masses are above the threshold of fermion 
pair excitations so that massive diquarks might eventually decay into two NG bosons. 3 The detailed numerical 
investigation of diquark masses as functions of the chemical potential is outside the scope of this paper and will be 
considered in a future publication. There, we are going to study also the masses of mesons in the environment of 
dense quark matter in the CFL phase. 
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Our earlier inve stig ations of the dispersion relations of diquarks in color superconducting quark matter, composed of u and d quarks 
(see the papers )15| ) , indicate that massive diquarks may occur as resonances which are heavily damped. However, the situation in the 
CFL phase might be different and needs a special consideration. 
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APPENDIX 



In the Nambu - Gorkov representation the inverse quark propagator matrix S 1 is given in ([8]) . Using the techniques, 
elaborated in [3, EH, EH, E3, EH > it is possible to obtain the following expressions for the matrix elements of the quark 

propagator matrix Sq = ( „ n ' c 12 

V '-'21, ^22 



Sn(x,y) 
Su(x,y) 
S2i(x,y) 
S22{x,y) 



d 4 q 



-iq(x-y) 



Q0-E+ - q + E 

-7 A+ + 



9o 2 - (^a) 2 



9o 2 - (E BA ) 



r7°A- , 



(2tt)4 



1 



:7 5 A- + 



1 



:2n) 



d V r -iq(x-y) 



(2tt)4 



- (E+ A f q 2 - (E BA ) 2 

qo-E 



7 5 A+ f , 



.?o 2 -(^a) 2 
Q0 + E+ 



7 5 A_ 



9o 2 - (^ 2 



where M 



3 erg, A^ 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



s 3 



M 2 and 



l± = t (l ± 7 (iq - lvl> ). Moreover, (£f A ) 2 = {E^ 2 + \A\ 2 B 2 , E ± = E ± fi, E = 

Z3 = X)a=2 5 7 t aAa- (I n these and other similar expressions, qo is a shorthand notation for qo + ie ■ sgn(go), where 
the limit e — > 0+ must be taken at the end of all calculations. This prescription correctly implements the role of [i as 
the chemical potential and preserves the causality of the theory.) 

It is clear from (|A.l|) - (IA.4j) that all color- and flavor dependences in the matrix elements Sn, S12, S21 and S22 arise 
only due to the matrix B. In the nine-dimensional space c x /, which is the direct production of color and flavor 
spaces, the two 9x9 matrices B and B 2 take the following forms 



B = 

















-1 











-1\ 




/ 2,0,0,0, 1,0,0,0, 1\ 













1 



















0,1,0,0,0,0,0,0,0 






















1 










0,0,1,0,0,0,0,0,0 







1 

























0,0,0,1,0,0,0,0,0 




-1 























-1 


, B 2 = 


1,0,0,0,2,0,0,0,1 

























1 







0,0,0,0,0,1,0,0,0 










1 






















0,0,0,0,0,0,1,0,0 



















1 













0,0,0,0,0,0,0,1,0 




-1 











-1 











J 




V 1,0,0,0,1,0,0,0,2/ 



(A.5) 



Thus, the quark propagator Sq may be thought as the matrix in the c X / space as well. Since it is an infinite series 
in powers of B, the diagonalization of the matrix B means the diagonalization of the propagator Sq in the c x f 
space. Let us define in the c x / space the following matrix (its rows are the ortonormal eigenvectors of the _B-matrix, 
corresponding to the eigenvalues p\ 



O = 





= P5 


= 1,P6 


= Pi 


= P8 




1,P9 


= -2): 




1 








- 


2 

V6 








7S 


\ 





1 

v/2 





1 

V2 






















1 

v/2 











1 






















1 

-J2 





* 1 




1 

V2 




















7f 







1 

V2 





1 

72 






















1 

V2 











1 

V2 






















-71 


* ? 




1 











1 








73 


/ 



(A.6) 



V 

Evidently, we have OO* = 1, detO=l and OB 2 O t =diag(l, 1, 1, 4) OBO t =diag(l, 1, 1, 1, 1, -1, -1, -1, -2). More- 
over, the matrix O diagonalizes the quark propagator Sq in the direct-product space c x /: OS^O 1 =diag(Si, Sq), 
where each of Si,...,Sg corresponds to a Nambu-Gorkov representation of a propagator of the quasiparticle with 
gap |A|, whereas S9 — to a quasiparticle with gap 2|A|. Hence, in the CLF phase all nine quasiparticles form an 
SU(3)-octet with gap |A| and an SU(3)-singlet with gap 2|A| (see also d, EU). 
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